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We present and discuss the results of a Monte-Carlo simulation of the phase transition in pure compact U(l) 
lattice gauge theory with Wilson action on a hypercubic lattice with periodic boundary conditions. The statistics 
are large enough to make a thorough analysis of the size dependence of the gap. In particular we find a non-zero 
latent heat in the infinite volume limit. We also find that the critical exponents v and a are consistent with the 
hyperscaling relation but confirm that the critical behavior is different from a conventional first-order transition. 



1. INTRODUCTION 

The interest about the nature of the phase tran- 
sition in compact 4D U(l) lattice gauge theory 
has been revived by the recent development of 
two new line of results. On the one hand, Kerler, 
Rebbi and Weber Q studied the critical proper- 
ties of the model by adding to the standard Wil- 
son action a coupling A controlling the density of 
monopoles. They concluded to the existence of 
a non-gaussian second-order critical point in the 
(/3, A) plane. Damm and Kerler Q are also in- 
vestigating whether the critical exponent of this 
transition is universal or changes with A. 

On the other hand, Jersak, Lang and Neuhaus 
m studied the compact U(l) gauge theory on lat- 
tices with sphere-like topology with a Wilson ac- 
tion extended by a coupling 7 of charge 2 



S' = /3 ^ cos Op -f 7 ^ cos 29^ 



p 



(1) 



They found no gap on these lattices for 7 < 0. 
They also made a thorough finite-size size scaling 
analysis of their data and concluded to the exis- 
tence, for 7 < 0, of a second-order transition with 
a non-gaussian continuum limit. 

Finally there is an investigation Q of the 
scaling behaviour of gauge-ball masses and of 
the static potential, which seems to confirm the 
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second-order nature of the transition also on lat- 
tices with periodic boundary conditions at 7 = 
-0.2 and 7 = -0.5. 

However we must note that the critical expo- 
nents of these two approaches are different. More- 
over there is always an apparent contradiction be- 
tween the simulations on lattices with spherical 
topology and on lattices with periodic boundary 
conditions since one observes a gap on the latter 
even when 7 < 0. One can fairly state that some 
confusion about the nature of the transition still 
persists. 

Therefore it is useful to reconsider the simula- 
tion of compact pure 4D U(l) lattice gauge the- 
ory with the standard Wilson action (7 = 0) and 
periodic boundary conditions provided that such 
a simulation fulfills two goals not met previously 

• give an estimation of the infinite volume 
limit of the gap which is observed on finite- 
size lattices with periodic boundary condi- 
tions. 

• make a careful modern finite-size scaling 
analysis of the bulk critical behavior (as has 
been done on lattices with sphere- like topol- 

ogy i). 

In order to estimate the infinite volume limit of 
the gap, one needs extrapolation formulas which 
require at least 3 parameters. Also the asymp- 
totic scaling formulas depend in general upon 3 
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parameters. Determining these parameters from 
the measurement of one observable only requires 
at least 6 to 7 data points. One needs more data 
points if one wants to have some chance of esti- 
mating the subleading correction terms. 

Therefore to be as systematic as possible the 
simulation has been done on 9 lattice sizes from 
L = 4 to L = 16. One could argue that these 
linear sizes are too small to reach the asymp- 
totic scaling regime. If this turns out to be in- 
deed the case, then these smaller lattice sizes are 
anyhow needed to determine the corrections to 
scaling which will be required even when data on 
larger lattices become available. 

2. SIMULATION 1 

Since accurate data were lacking for many of 
the above lattice sizes, a first simulation has been 
done to determine their pseudo-critical coupling 
to about 1 part in 10""*. In order to reach this 
accuracy a scanning of the pseudo-critical regions 
was done with a step A/? — 10^** and 10^ itera- 
tions at each coupling constant /3. 

We used standard histogramming techniques to 
locate the double-peak structure found in these 
simulations which is usually characteristic of a 
first-order transition. A pseudo-critical coupling 
was defined as the coupling for which peaks have 
equal statistical weight. The latent heat was de- 
fined as the gap between the two peaks at this 
pseudo-critical coupling. 

The details of this simulation are described in 
1^ and we only give here a short summary of the 
results. 

A three-paramater fit to the gap Ae{L) of the 
form 

Ae{L) ^ Ae{oo) + aL-'' (2) 

gives Ae(oo) = 0.014(5) and b = 1.03(16). We 
note that the exponent b is not consistent with 
the first-order prediction (at least in Potts mod- 
els) 6 = D = 4|||,|l^. Fixing 6 to 1 reproduces 
the data quite well but with an infinite volume 
limit of the latent heat which is definitely differ- 
ent from zero. However it should be stressed that 
it is very difficult to constrain the functional form 
of a three-parameter fit to the gap data. Indeed 



an exponential fit 

Ae(L) = Ae(oo) + a exp(-6i) (3) 

reproduces the data as well with Ae(oo) = 
0.0278(15). 

An asymptotic finite-size scaling analysis of the 
pseudo-critical couplings 

(3ciL) = f3ci(x,)+aL-i (4) 

yields the results bc{oo) = 1.01132(10) and v = 
0.326(8). 

We have also checked the scaling of the maxima 
of the specific heat for lattice sizes in the range 
i = 4 — 12. A two-parameter ansatz 

Cv,max{L) = aL~^^ (5) 

gives the independent determination i> — 0.330(2) 
but with a rather high ~ 3 which hints at 
the need of correction terms to the asymptotic 
formula. 

All these results taken together confirm the 
rather paradoxical nature of the U(l) phase tran- 
sition. The critical exponents are completely con- 
sistent with a second-order phase transition but 
with an index v « 0.33 which is different from the 
value, v « 0.36, quoted in This discrepancy 

raises the suspicion about universality at different 
values of 7. On the other hand any reasonable fit- 
ting ansatz to the gap data yields a non-zero value 
of the latent heat in the infinite volume limit. But 
again the approach to this limit is different from 
the asymptotic formula expected within the de- 
scription of first-order transitions in the double 
gaussian approximation This disagreement 

might mean that the asymptotic regime is not yet 
reached with lattice sizes up to L = 16. 

3. SIMULATION 2 

Going to larger lattices is impracticable with 
local algorithms since already we could not over- 
come the hysteresis on the 16"* lattice with 10^ 
iterations. However it is possible to attack the 
problem indirectly by increasing the statistics on 
the smaller lattices so as to make a full finite 
size analysis including corrections to scaling. The 
comparison of the finite size scaling of several cu- 
mulants can unravel the systematic errors in the 
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critical exponents induced by the corrections to 
scaling. Combined fits can reveal whether the 
critical exponents vary when excluding the small- 
est lattice sizes. 

Therefore we have made a second simulation at 
3 to 5 coupling constants selected in the pseudo- 
critical interval determined in simulation 1 at 
each lattice size, except L = 16. 10^ iterations 
have been done at each coupling constant. These 
10® iterations were divided in two independent 
runs, 5 X 10^ sweeps each, respectively from a 
hot start and a cold start, using different random 
generators. 

The data analysis, which is not yet completed, 
makes an extensive use of the reweighting tech- 
nique | pl| . All independent runs on the same lat- 
tice size are used as independent samples at the 
same f3. The total amount of statistics that we get 
is quite comparable to many of the Monte-Carlo 
simulations of 3D spin models. 

We are making a finite-size scaling analysis 
of 3 cumulants: the specific heat per plaquette 
C„, the Binder Cumulant U4 and the second cu- 
mulant U2- We are also adding an analysis of 
their derivatives dCv/d(3, dU4/d(3, dU2/df3. These 
6 cumulants are algebraically independent and 
the position of their extrema defines a pseudo- 
critical coupling if located in the scaling region. 
The value of each pseudo-critical coupling is de- 
termined independently for every run by mini- 
mizing the corresponding reweighted histogram 
with respect to (3. Finally we take the statisti- 
cal average over all runs at each lattice size. The 
caveat of the method, and the limiting factor of 
its applicability, is to ensure that the Ferrenberg- 
Swendsen technique remains valid throughout the 
minimization process. 



Preliminary results from independent asymp- 
totic finite size scaling fits to the various defini- 
tions of the pseudo-critical couplings show devia- 
tions « 2% among the values of the critical expo- 
nent J/. Even if these deviations are much larger 
than the statistical errors, it will be very difficult 
to extract the corrections to scaling. 
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4. CONCLUSION 

Preliminary results of the second simulation 
confirm the non-zero value of the latent heat in 
the infinite volume limit. Fits for Ae(cx)) with 
Eq. H and Eq. || are completely consistent between 
both simulations. However the value of the pa- 
rameter b in Eq. ^ increases to « 1.40 and the 
corresponding is much larger than for the ex- 
ponential fit which becomes highly favored. 



